Phenomenology of a semi-Dirac semi-Weyl semi-metal 



S. Bancrjec^ and W. E. Pickett^ 
^Department of Physics, University of Califor-nia, Davis, CA 95616 
(Dated: July 25, 2012) 

We extend the study of fermionic particle-hole symmetric semi-Dirac (alternatively, semi-Weyo) 
dispersion of quasiparticles, ek = ±\/{k'^/2m)'^ + (vky)^) = zL£o^yK* + K'^ in dimensionless units, 
discovered computationally in oxide heterostructures by Pardo and collaborators. This unique sys- 
tem a highly anisotropic sister phase of both (symmetric) graphene and what has become known as 
a Weyl semimetal, with < >^/^« v independent of energy, and < >^''^oc nC^^"^ ^ being very 
strongly dependent on energy (e) and depending only on the effective mass m. Each of these systems 
is distinguished by bands touching (alternatively, crossing) at a point Fermi surface, with one con- 
sequence being that for this semi-Dirac system the ratio IXor-t/Xspl of orbital to spin susceptibilities 
diverges at low doping. We extend the study of the low-energy behavior of the semi-Dirac system, 
finding the plasmon frequency to be highly anisotropic while the Hall coefficient scales with carrier 
density in the usual manner. The Faraday rotation behavior is also reported. For Klein tunneling for 
normal incidence on an arbitrarily oriented barrier, the kinetic energy mixes both linear (massless) 
and quadratic (massive) contributions depending on orientation. Analogous to graphene, perfect 
transmission occurs under resonant conditions, except for the specific orientation that eliminates 
massless dispersion. Comparisons of the semi-Dirac system are made throughout with both other 
types of point Fermi surface systems. 

PACS numbers: 



I. INTRODUCTION 



The isolation of single layers of graphite 
(graphene) with its unique linear (massless Dirac, 
properly ealled Weyl) low energy band structure has 
become, within only a few years, a heavily studied 
phenomenouiii^ The appearance of unanticipated 
new features in band structures, which generally 
have far-reaching implications, have in the past in- 
cluded half metallic ferromagnets and compensated 
half metals ("half metallic antiferromagnets" ) , and 
more recently topological insulators j^ii Each of these 
systems provide the promise of not only new phys- 
ical phenomena but also new applications of their 
unconventional properties. 

Another key feature of graphene is the point Fermi 
surface aspect. The touching (or crossing) of bands 
is accompanied by a gap throughout the rest of the 
Brillouin zone that pins the Fermi level (E^?) in 
the intrinsic material to lie precisely at the point 
of crossing - the point Fermi surface (two of them in 
graphene) . This point Fermi surface aspect has been 
well studied^ in conventional zero gap semiconduc- 
tors where a touching of the valence band maximum 
and conduction band minimum is symmetry deter- 
mined and occurs at a high symmetry point. The di- 
electric susceptibility of such a system is anomalous^ 

- neither metallic nor semiconducting in character 

— and unusual consequences of the touching bands 



and residual Coulomb interaction promise unusual 
phases, such as excitonic condensates including ex- 
citonic superconductors and excitonic insulators. 

The linear dispersion at the zone boundary in 
graphene has been known for many decades; it took 
the ability to prepare the delicate material and per- 
form a variety of experiments to ignite interest. 
There are quasilinear (and potentially truly linear) 
band structure features in certain materials, viz. 
skutteruditesji that have been known for some time 
and with recent developments^ may attract new at- 
tention. To actually discover a feature in a band 
structure that provides the quasiparticle dispersion 
of a new and unexpected type is rare, and the discov- 
ery of a semi-Dirac dispersion pinned to the Fermi 
energy is a very recent example. 

Pardo and one of the author o^'^° reported such 
a finding in ultrathin (001) VO2 layers embedded in 
Ti02. This new point Fermi surface system, dubbed 
'semi-Dirac,' is a hybrid of conventional and uncon- 
ventional: dispersion is linear ("massless", Dirac- 
Weyl) in one of the directions of the two-dimensional 
(2D) layer, and is conventional quadratic ("mas- 
sive") in the perpendicular direction. At directions 
between the axes the dispersion is intermediate and 
highly direction-dependent. Interest in this unique, 
maximally anisotropic, dispersion arises for several 
reasons. The (topologically determined pinning at 
the) point Fermi surface is itself of interest. The 
highly anisotropic dispersion (from massive to mass- 



less depending on angle) is unique to this system. 
The fact that it arises in an oxide nanostructure of 
the general type that is grown and studied regularly 
these days also strengthens the promise of appli- 
cations. Another layered superstructure, a double 
cell layer of Ti3SiC2 embedded in SiC, has displayed 
a point Fermi surface, but the dispersion is of the 
convention typcj^ As alluded to above, an unusual 
point Fermi surface at zero momentum, with linear 
bands degenerate with quadratic bands, has been 
discovered in the skutterudite class of semimetals4 

Such a spectrum had been noted earlier in differ- 
ent contexts. Volovik obtained such a spectrum at 
the point node in the A-phase of superfluid '^He [l^l 
and studied its topological robustness>i^ More rele- 
vant to solids was the discovery by Montambaux's 
group of this spectrum in a graphene-like model.— 
The model has a broken symmetry such that hop- 
ping to two nearest neighbors is t but to the third 
neighbor is t' . When t' differs from t, the graphene 
"Dirac points" wander away from the K and K' 
points, and at t' = 2t they merge, resulting in the 
semi-Dirac spectrum. This group began a study of 
low energy properties of such a system^i^ which was 
continued by Bancrjee et alr^ and will be extended 
in the present paper. 

In this paper we first provide results for the Hall 
coefficient and plasma frequency versus doping level, 
finding some new behavior along with some some- 
what conventional results. In the final section we 
provide selected results for Klein tunneling of semi- 
Dirac particles, a problem that acquires extra rich- 
ness due to the variable angle of the barrier with 
respect to the anisotropic dispersion. 

II. SEMI-DIRAC DISPERSION 

SemiDirac dispersion is quadratic along one sym- 
metry direction in the Brillouin zone and linear along 
the direction perpendicular to it: massless Dirac(i.e. 
Weyl). Choosing and ky to be the momentum 
variables and taking h=l except occasionally for 
clarity, the semi-Dirac dispersion is given by: 




where the effective mass m applies along k^ and v is 
the velocity along A:j,(the massless direction). For in- 
termediate angles /3 = a.icta.n{ky / ky) , the dispersion 
is of an entirely new type. Two natural scales are 
introduced, one for the momentum and the other 



for the energy: p = 2mv (momentum scale) and 

2 

^0 ~ 2m ~ '^P'^- (Untidy factors of 2 appear be- 
cause of the clash between the natural classical ^pv 
and relativistic pv units for energy.) One can then 

define the dimensionless momenta Ky = and 

^ p 

hk 

Ky — in terms of which the semi-Dirac disper- 
sion given by Eq. [1] becomes 

ek = ±£o/^JT^. (2) 

The corresponding velocity = Vk£k can be scaled 
to a dimensionless form Vk using 

Vk = — = (3) 

V 

Figure [T] shows semi-Dirac Fermi surfaces as well as 
contour plots of Vk- 

We first compute (t;^) and (w^), which arc the av- 
erages of the Fermi surface velocity vp = ((w^) + 
(Vy))^ for the semi-Dirac dispersion in the non- 
relativistic and the relativistic directions respec- 
tively, which will prove to be useful later, and will 
also give the semiclassical conductivity tensor aap = 
e'^TD{e){vaVp) . They are defined as follows 

{vavpie)) = Vavi36{ek - e)/ 5{£k - e) (4) 

k k 

27r^D{e)J ' \vk\ ' 

where D{e) is the density of states. For semi- 
Dirac dispersion the density of states was obtained 
earlier— as 

V2me 2m [T 
D{e)=Ii — 5 — / — , (5) 

with proportionality coefficient mjv^. The inte- 
gral I\ is given by 

h= I ds(l - s")"^ « 1.3110. (6) 

Bacsi et al. have studied the quantum critical ex- 
ponents of point Fermi surface semimetals^^ with 
D{e) (X \s\^ for a continuous range of r including 
this r = 1/2 case. 

The squared Fermi velocities for scmiDirac disper- 
sion arc obtained as 

{vD ^'f-^^ 1.3.^ (7b) 
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Figure b 




FIG. 1: Fig. a: Fermi-surfaces of semi-Dirac dispersion, 
along with arrows representing Vk- The length of an ar- 
row is proportional to the magnitude of Vk- As can be 
seen from the figure, the arrow-length is constant along 
the Ky axis indicating a constant velocity in the rela- 
tivistic (j/) direction. Vks are all normal to the constant 
Fermi energy contours, as they should be. Fig. b: The 
surface and contour plot of the magnitude of Fk- The 
magnitude is constant in the y direction as opposed to 
the monotonically changing values in the non-relativistic 
(x) direction, with rapid variation of other directions of 
propagation. 



note that the former involves only m, the latter only 
V. The integrals I2 and /s are given by 



Jo 

Jo (l-.s4)i 



0.8740, (8a) 
i 0.3595, (8b) 



Thus the ratio of (v^) to (Vy) scales as e/eo, which 
reflects the extreme anisotropy at small doping. For 
the VO2 system where semi-Dirac dispersion was 



discovered only very small doping levels will re- 
main within the energy range represented by the 
semi-Dirac dispersion (^ ~ 10~^) but we consider 
more general cases. 



III. FARADAY ROTATION IN THE 
CONTEXT OF THE SEMI-DIRAC SYSTEM 

A. The semiclassical equation of motion 

The behavior of point Fermi surface semimetals 
in a magnetic field has stimulated lively interest due 
to unusual quantum Hall effect behavior, with the 
case of graphene having been reviewed recently by 
GocrbigJ^ The semiclassical equation of motion of 
an electron in a magnetic field B is given by 



n— 

dt 



^-vl X B. 
c 



(9) 



Using Eq. Q for vl in Eq. ([9]), one obtains the fol- 
lowing expressions 



dK^ 

dKy 

dt 



= -UJoKy, 



(10a) 
(10b) 



where K,j: and Ky are the dimensionless variables 
associated with momentum introduced before, and 
loq is given by 



Wo 



eBv eB Eq 
cs mc e ' 



(11) 



where B is the magnetic field, and e, the Fermi en- 
ergy. The Fermi surface orbiting frequency diverges 
as the doping level decreases; the Fermi surface or- 
bit length goes smoothly to zero whereas the mean 
velocity remains finite. Eliminating Ky from Eqs. 
()10a|) and (jlOb|) . the following differential equation 
is obtained 



(PK^ 
dt^ 



(12) 



In order to solve this second order differential equa- 
tion, we multiply both sides of the equation by Kx 
(Kx denotes the time derivative of Kx). Both the 
right and the left sides of the equation can then be 
written as a total derivatives, which can be inte- 
grated to give 



K^. 



(13) 



3 



where the constant C can be determined from 
the condition that = when = -ftT^j^max- 
-^a.max — {sjeof-l'^ corresponds to Ky = 0, from 
the semiDirac dispersion given by Eq. [51 and the 
rest follows from Eq. (|10ap . Hence Eq. [T3] becomes 



C. Faraday Rotation 



The Faraday rotation angle is given by the 
expressions^ 



a:, max x 



(14) 



Integrating the above equation (numerically) one 
can get K^, as a function of time. Once is known, 
Ky can be obtained from Eq. llObI The differential 
equation for the cyclotron orbit is obtained by divid- 
ing Eq. llObl bv Eq. llOal Solving for that, we obtain 
the semi-Dirac constant energy contour as an expres- 
sion for the cyclotron orbit, which is expected, since 
the energy of an electron does not change when it 
moves under the influence of magnetic field. 



B. The cyclotron frequency 



Eq. [T3]can be integrated using the limit — i^x.max 
to A'a;.,„ax for the variable to obtain the time 
period. The result for the time period (T) thus ob- 
tained is 



4/i 



= 4/i Ve/eo , 



(15) 



where 1\ is given by Eq. [51 From Eq. 1151 the fun- 
damental semi-Dirac cyclotron frequency il, 
obtained as 



= f IS 



(16) 



The cyclotron frequencies for the parabolic and 
the linear dispersion cases are given by (^%^ = .2^) 



and 



respectively (/is is the Bohr magneton). 



Comparing with Eq. [12] we see that the cyclotron 
frequencies for all the three cases(the parabolic, lin- 
ear, and semi-Dirac) depend linearly on the mag- 
netic field. The cyclotron frequency is independent 
of the Fermi energy for parabolic dispersion, whereas 
it varies as e~2 for the semiDirac dispersion and as 
for the linear Dirac dispersion. One important 
aspect of the semi-Dirac dispersion is that the semi- 
Dirac dispersion being anisotropic in the momentum 
space can have harmonics of the fundamental cy- 
clotron frequency given by Eq. 1161 This feature is 
absent in the Dirac or the two dimensional parabolic 
dispersion where the energy momentum dispersion 
is isotropic giving rise to only one value for the cy- 
clotron frequency. 



B(uj,B) = Zo/,(w)i?e[t7,.y(w,B)], (17) 

where is the impedance of the vacuum, /g is the 
spectrally featureless function specific to the sub- 
strate, and Oxy is the dynamic Hall conductivity. 
According to the Drude formula the dynamic Hall 
conductivity is given byi^ 



-21? 



'xy 



(18) 



where 2? is the Drude weight, given by P = 
|e2£)(e)(u2). Taking the real part of Eq. [Hj and 
using it in Eq. [T7| we obtain 



(w,B) = /(w), (19) 



where /(w) is given by 



/(c) 



(20) 



Extremizing /(w) and inserting the resulting ex- 
pression for /(w) in Eq. [T9| wc obtain the following 
expression for the maximum value of the Faraday 
rotation angle 



(w,5) 



-Zofs{uj)'Dujc 



I) 



2)' 



(21) 



The Drude weight T) ^ e for Dirac disper- 
sion(since D{£) ^ e, and (w^) is a constant). The 
Dirac cyclotron frequency ijJc^ ■ Hence the prod- 
uct "Dujc that appears in the numerator of Eq. [21] is 
independent of the doping level for Dirac dispersion. 
For semi-Dirac dispersion, T) ^ , which follows 
from the fact that the product D{e){v^) ~ D{e){Vy), 
where Vy is the speed in the relativistic direction, and 



that D{e){i 



£ 2 . The last step follows by com- 



bining Eq. [5] and Eq. I7bl For the same dispersion 
ujc ^ e (From Eq. [T6)) . Hence, like Dirac dis- 
persion, TDuJc for the semi-Dirac dispersion is inde- 
pendent of the doping energy. For two dimensional 
parabolic dispersion, ujc is independent of the dop- 
ing energy, but T) £. Hence "DuJc depends on the 
doping energy. This is a significant difference when 
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compared to the Dirac and the semi-Dirac disper- 
sion. 

For Dirac and semi-Dirac systems the dependence 
of the Faraday angle on the doping level arises from 
the term uJcT in the denominator of Eq. 1211 whereas 
the numerator is independent of doping. For those 
dispersions one can fine tune the Fermi energy to 
obtain a large value of the Faraday angle by bring- 
ing the term lOcT close to three, so that the term 
(w^T^ -1-1)2 — 2 appearing in denominator goes to 
zero causing a significant value for the Faraday an- 
gle. 



IV. HALL COEFFICIENT 

According to semiclassical Bloch-Boltzmann 
transport theory, the Hall coefficient of a two 
dimensional Fermi liquid (in the x — y plane) is^S 



R 



H 



R 



H 
xyz 



S]i.^;.(k)[v(k) X V(k)]3z;,(k)(^ ), 



(22) 



Due to the algebraic complexity of the first and sec- 
ond derivatives of , this expression is formally un- 
wieldy. We show however that general properties of 
this expression lead to a simple and familiar result 
for Rh- 

The numerator of Eq. is the area Ay spanned 
by the velocity vector over the Fermi surface^!. In 
the zero temperature limit each term in the denom- 
inator reduces to a line integral along the Fermi- 
surface. The carrier density n is proportional to the 
area swept by the vector k over the Fermi surface, 
which is the area Aps enclosed by the Fermi surface. 
Hence the quantity R^n is given by: 



Ay A 



FS 



J <■ Vk J ' Ufc 



(23) 



Using the fact that the gradient VkS is perpendic- 
ular to the vector line element dki along the Fermi 
surface, so that the dot product between them is 
zero, the denominator of Eq. reduces to 



(24) 



y dki — y dki ~ ~ dkyVx <j) dk^Vy 

Using Eq. ^ in Eq. ^ we obtain 
AyAps 



R.^n 



dkyVx j> dkxVy 



(25) 



R^n as given by Eq. is unity for the semi- 
Dirac dispersion. This result can be argued directly 
from equation Eq. (|25p in the following way. The 
semi-Dirac dispersion is symmetric both in the x and 
the y directions. Hence we can restrict the limits 
of the integrals appearing in Eq. (^5]) to the first 
quadrant. For the first term in the denominator of 
Eq. (psj) . carrying out the integration by parts one 
obtains: 



dky 



-kyVx\l 



f 



dvx ky 



(26) 



i and / correspond to the points on the Fermi surface 
with ky — and k^ ~ respectively. The boundary 
terms in Eq. (|26p at i and / are zero because ky and 
the X component of the gradient at the semi-Dirac 
Fermi surface vanish at i and / respectively. Using 
the above reasoning the first term in the denomi- 
nator of Eq. (|^ is changed to / dvxky. Making 
use of this along with the definition of area under a 
curve(for the terms in the numerator), Eq. (j25p can 
be written as 



R"n 



J dkxky J dvxVy 
J dvxky J dkxVy 



(27) 



Vy for the semiDirac dispersion evaluated on the 
Fermi surface turns out to be proportional to ky as 
can be seen from Eq. ©. Hence it is observed that 
in Eq. (|27p the numerator and the denominator are 
equal except for a minus sign. That explains why 
we obtain R^n — —1 for the semiDirac dispersion. 
Incidentally, Vy is proportional to ky for the Dirac 
and the parabolic dispersion relations. Hence, R^n 
is equal to —1 for those dispersions too. So it can be 
said that the Hall coefficient times the carrier den- 
sity is a topologically invariant quantity for a certain 
class of band structures, reminiscent of the geomet- 
rical representation of Ongi^ 



V. PLASMON FREQUENCY 

The plasmon frequency for the semiDirac system 
can be computed by setting the random phase ap- 
proximation expression for the dielectric constant 

e(q,c<.) = l-i'(q)Xo(q,w) (28) 

to zeroi22i2^ xoi'i,'^) is the polarizability and v(q) 
is the Fourier transform of the Coulomb potential. 
Xo(q, '^) is given by the Lindhard expression 

d'k /(£k) - /(ek+q) 



xo(q,t^) 



(27r)2 



(29) 
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Expanding ^k+q in Eq.[2Hlfor small q (we treat only 
this regime), the numerator in Eq. [21] takes the fol- 
lowing form at low temperature 

/(^k) - /(£k+q) = Vk ■ qS{ei, - e). (30) 

Expanding the denominator as well, Eg. 1291 becomes 

Xo q,w = / — 1 + (5 £k-e . 31 

The Coulomb potential w(q) in two dimensions is 

«(q) , (32) 



where q ~ y Ix + 1y ■• ^nd k is the background di- 
electric constant of the medium. Using Eq. [21] and 
Eq. [32]in Eq. and setting e{q,uj) = 0, the plas- 
mon frequency is 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

angle between the plasmon wave vector q and tlie non relativistic direction 



FIG. 2: Angular dependence of the function F 



8/3 e^geo 



(33) 



where F(0) is given by 



F{0) = ^^cosH + lr^^sinH), (34) 

and I2, I3 are given by Eg . [8al and in Eq. [8b] respec- 
tively. Eq is the energy scale defined earlier. 6 de- 
notes the angle that the plasmon wave- vector makes 
with the non-relativistic axis of the semi-Dirac 
dispersion. Recall that the Fermi energy variable is 
defined as ^ = Up cc ^/q is characteristic of a 
two-dimensional system. 

The function F{6) is plotted against 9 in Fig. [2] 
Using Eq. [5] for the semiDirac density of states and 
Eq. [7i and Eq. [7b] for the mean square Fermi veloc- 
ities. Eg. [551 reduces to 



e^qhD{e 



-{{vl)cos'e+{vl)sinH){i^) 



The plasmon frequency is highly anisotropic and 
reaches its maximum along the relativistic direc- 
tion, which could be a signature characteristic of s 
semiDirac system. 



VI. MAGNETIC SUSCEPTIBILITY 

In this section we consider the magnetic suscepti- 
bilities for the semi-Dirac dispersion. The Pauli spin 



susceptibility is given by 



(36) 



where D(e) is the density of states. Using Eq. [5]for 
the semi-Dirac density of states Eq. [36] reduces to 



(37) 



where ^ is the same dimensionless variable related 
to the Fermi energy appearing in the previous sec- 
tion. For a non-interacting Fermi liquid the orbital 
susceptibility is given by2i 



Xorb/ fJ^B 



(9^gk 5^£k 
dk"^ dky 



3 0£k 9^£k 9gk 5^£k ^1.. _ s 

2 ^ dk, dk^dkl dky dkydkl 

Using Eq. [2] for £k in Eq. [38] and doing the integral 
we obtain 



Xorb/ ^J-^ 



2V2L 



4 m^v 



37r3 £5 

where the integral I4 is given by 

, -33ai" + 41a*^ - 

14= da ■, . 

Jo (l-a4)^ 



(39) 



(40) 



Evaluating the numerical value for I4 and using the 
dimensionless variable ^, Eq. [39] reduces to 



Xorbl IJ.'b - - 



0.0798m jeo 



(41) 
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We observe that the orbital susceptibihty for the 
semi-Dirac system is always diamagnetic. The ab- 
solute value of the ratio of the spin to the orbital 
susceptibilities (the ratio of Eq. [37]to Eq. (HTI) ') of 
the semi-Dirac dispersion is given by 



Xorb 



looe 



(42) 



Hence orbital magnetic susceptibility for the semi- 
Dirac dispersion dominates the spin susceptibility at 
low energy. This result is distinct qualitatively from 
both the Dirac and the parabolic dispersion cases. 
For the doped Dirac dispersion the orbital suscep- 
tibility vanishes identically. For conventional two 
dimensional parabolic dispersion the orbital suscep- 
tibility is calculated using Eq. I38[ and turns out to 
be 6tt times smaller than its paramagnetic suscep- 
tibility. Hence the unusually large orbital suscepti- 
bility can be considered a distinctive feature of the 
semi-Dirac dispersion. 



VII. HEAT CAPACITY 



electron gas is first outlined in the following. The 
parabolic three dimensional Hamiltonian is given by 
^fparaboiic = ^{pI + pI + pI), SO it foUows that 
agpa^aboiic ^ 2± [where i = x,y,z]. Hence TJparaboiic 
can be written as 



parabolic 



1 / ^-^parabolic 



(44) 



^^parabolic , ^-^parabolic \ 

+Pv %I ^ )• 



dpy 



dpz 



By the equipartition theorem, the ensemble average 
ofeachofp.^^i^,p,^%^,andp,^^i= 

is ksT^ Hence taking the ensemble average of the 
Hamiltonian in Eq. |44l one obtains 



< ^^parabolic >— I^^bT. 



(45) 



The derivative of < iJparaboiic > with respect to T 
gives the heat capacity as 

Next, the classical semi-Dirac Hamiltonian is 
given by 



We show here how the heat capacity for the 
non-interacting two-dimensional semi-Dirac electron 
gas is similar to that of the t/iree-dimensional non- 
interacting electron gas with the parabolic energy- 
momentum dispersion at both the low and the high 
temperature ends. The similarity becomes equality 
at high temperature. Relative to the natural energy 
scale Eq introduced at the beginning, the low and the 
high temperatures can be considered. The low tem- 
perature heat capacity per particle for the semiDirac 
dispersion is : 



2/i 



mk\,T 



(43) 



which is calculated using Sommerfeld expansion=i(/i 
is given in Eq. It is observed that the heat ca- 
pacity in Eq. |43] is proportional to D{e) oc ^/e, as 
it must be because depends only on the spec- 
trum of energy levels. A similar type of dependence 
with energy is observed for the three dimensional 
electron gas with the parabolic energy-momentum 
dispersion. The difference between them is in the 
prefactors. This difference disappears quite nicely 
in the high temperature end as is shown in the 
following. At high temperature, the heat capac- 
ity for the three dimensional electron gas is given 
by ffefi. In order to emphasize a technique that 
will be used for the semi-Dirac problem, a deriva- 
tion of the above result for the three-dimensional 



Am? 



(46) 



Taking the derivatives of Hgu with respect to Px and 
Py gives, in spite of its complex form, the analogous 
expression 



1 dH,u 
2^^ dpx 



Py- 



dHsD 
dpy 



(47) 



In the same way as before, by the equipartition theo- 
rem, the averages of each of Px^§^,Py^§^ is fcsT. 
Hence the ensemble average of HsD is given by 

< i/,D \kBT + keT = hgT (48) 

thus Cy ~ ^ks for semi-Dirac dispersion in the high 
T limit. This result is exactly that of a three dimen- 
sional non-interacting gas with parabolic dispersion. 

This rather unexpected result can also be obtained 
directly starting from the Boltzmann distribution. 
In the low temperature limit the semi-Dirac heat 
capacity has the same T dependence as the non- 
interacting three dimensional parabolic system. In 
the high temperature end of the spectrum the heat 
capacities are identical. Hence a two dimensional 
semi-Dirac system effectively behaves as a three di- 
mensional system so far as heat capacities are con- 
cerned. The appearance of this third degree of free- 
dom can have potential technological applications. 
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For example, a semi-Dirac nanostructure could be 
used as an efficient heat sink. More generally, a semi- 
Dirac system can function quite differently com- 
pared to other two dimensional systems for thermal 
management as well as for many other applications. 

VIII. KLEIN PARADOX 

The Klein paradox is the name given to the phe- 
nomenon of the complete transmission of a parti- 
cle at selected energies or geometric configurations 
through a potential barrier even when the barrier 
is arbitrarily high. For the conventional tunnel- 
ing problem, the probability of transmission de- 
creases exponentially with the height and thick- 
ness of the barrier. In order for Klein tunneling to 
take place, there must be hole states having nega- 
tive energies available to promote tunneling. The 
positive potential in the barrier region raises the 
hole states, making them available. For 'relativis- 
tic' Dirac-Weyl dispersion (as in graphene) Katsnel- 
son and collaborators^! have shown that Klein tun- 
neling can occur and that transmission is unusually 
robust at near-normal incidence. Klein tunneling 
is also possible in conventional (massive) zero-gap 
semiconductors including double-layer graphene?^ 
with an angular behavior that is distinct from that of 
graphene. Klein tunneling therefore is expected for 
particles with semiDirac dispersion, but there should 
be many distinctions. The low-energy Hamiltonian 
corresponding to the semiDirac dispersion can be 
taken asi^ 

H = vpyTs + j^Ti, (49) 

where the r's are the Pauli matrices in orbital space 
and Px{y) are the momentum operators. 

A. Rotation of the Frame 

The semi-Dirac system is (highly) anisotropic. 
The potential barrier can be oriented at an arbi- 
trary angle with respect to the x, y axes, after which 
one might consider a particle impinging on the bar- 
rier from another arbitrary angle. This extension 
from the isotropic systems of graphene or zero-gap 
semiconductors leads to a rather complicated tun- 
neling problem that could form the basis of a sep- 
arate study. To keep the algebra and the physi- 
cal picture as simple as possible, we consider only 
the special case of normal incidence of a semi-Dirac 



y 




FIG. 3: The top view of the potential barrier is shown. 
It extends infinitely in one direction (7) direction), but 
limited to a spatial length d in the orthogonal direc- 
tion (^ direction), which makes an angle a with the non- 
relativistic direction. An electron with energy E is inci- 
dent normally on the potential, i.e, along the ^ direction. 



quasi-particle onto a potential barrier of width d, 
which is inclined at an angle ^ + a with respect 
to the a:;(nonrelativistic) axis as shown in FigjSl A 
set of orthogonal axes ^ and rj with respect to the 
barrier arc defined, f) is the direction along which 
the potential is infinitely extended. The electron is 
incident on the potential along ^, which makes an 
angle a with respect to the x axis. The barrier has 
thickness d along the ^ axis. We work in the regime 
where the energy of the incident semi-Dirac particle 
is small compared to the barrier potential. There 
are three real space regions: to the left of the bar- 
rier where the potential is zero; within the barrier 
with positive potential V; and to the right of the 
barrier where again the potential vanishes. We refer 
to these regions as /, //, and ///, respectively, and 
the wavefunctions are denoted by ^'/, v]>jj, "^m, re- 
spectively. The momentum operators along the x 
and the y (rclativistic) directions can be written in 
terms of the variables ^ and 77 as follows: 

Px — p^cosa ~ PrjSma (50) 
Py ~ sina cosa, 

where are the corresponding momentum oper- 
ators given by —id/d^{'q). Since we are considering 
incidence normal to the barrier, it is straightforward 
to show that the ?y degree of freedom can be elimi- 
nated from the problem. The Hamiltonian in Eq. 031 
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takes the following form: 

p? 2 

H = wpe sin araH ^ cos ari (51) 

2m 

This transformed kinetic "Hamiltonian" has both 
linear (massless) and quadratic (massive) contribu- 
tions, governed an increased mass m^ ~ m/cos^a 
and a decreased velocity Va = vsina. Thus the 
orientation of the barrier allows the tuning of the 
relative amounts of linear and quadratic dispersion. 
In the limits a = and 7r/2, the problem reverts 
to the problem for zero-gap semiconductors and for 
graphene, respectively. 

For a value of a between these limits the forward 
propagating wave, which is of the form e'*^^ times a 
spinor, is still an admissible eigenstate of the Hamil- 
tonian. Operating on the planewave with the Hamil- 
tonian in Eq. [21] gives an expression that can be 
written as 



Hk — ufcsina[r3 -|- ri tan0]. 



where 



tan^ 



cos^ a k 



k 



sin a 2mv 2toq,Uq, pa 



(52) 



(53) 



Thus tanO reflects the magnitude of the particle mo- 
mentum relative to the scaled semi-Dirac momen- 
tum Pa = 2inaVa- Whcu k goes to — k as is the case 
when one considers the backward propagating wave 
e"**^^, aside from the positive multiplicative factor 
vk which changes sign, the Hamiltonian in Eq. ISTl 
changes from + ti tan 9 to — [ra — ri tan 6] . The 
corresponding eigensystems are given for quick ref- 
erence in the Appendix. 



B. Derivation of the Resonance Condition 

The time independent Schrodinger equation in a 
given potential can be written as 



hi' = {E- F)V', 



(54) 



where h is the part of the Hamiltonian without the 
potential V. In regions / and /// {E — V) is pos- 
itive, and the positive eigenvalue form of the solu- 
tion as given by Eq. IA2al in the Appendix for the 
forward propagating wave and by Eq. IA4al for the 
backward propagating wave need to be considered 



in those regions. In region //, V being much larger 
than E results in {E — V) being negative. Hence the 
negative eigenvalue solutions as given by Eq. IA2bl 
and Eq. IA4bl appearing in the appendix are of im- 
portance in that region. Momenta in regions / and 
/// are equal, denoted by fci, and denoted by ^2 in 
region //. ki and ^2 are given by 



vki sinQ!(cos0i) 
vk2 sinQ!(cos 62) 



E, 

V-E, 



where 0i and O2 are given by 



tan0i(2 



cos^ a ki 



(2) 



fcl 



(2) 



(55a) 
(55b) 



(56) 



sin a 2mv pa 
Finally, the wave functions in the three regions are 



_ fc,^ /cos(0i/2) 
' " ^ Uin(ei/2) 



(57) 



, ( sin(0i/2)\ _^ < ^ < 

^ _ ^ psfec f sin(6'2/2) 
^" - (,-cos(^?2/2) 

*--^-'^^^KS(V2j)'^<^<- 

where r,ti,ri and ^2 are constants determined by 
matching. The absolute square of t2 gives the trans- 
mission coefficient. Matching the wave functions at 
the boundaries y = and y = d, one obtains for the 
transmission 



t2' = 



(sin 6*2 cos 02 cos 6*1 ) 



A2 +B2 _ 2 AB cos k2d' 
where A and B are given by: 



(58) 



A = [sin((6i2 -6ii)/2)cos6i2 (59) 

+ sin(6'2 + 6'i)/2] cos((6'2 - 0i)/2) 
B = sin02sin2((02+ei)/2) 



It can be shown that when 

cos k2d = 1 



(60) 



the denominator in Eq. [58] becomes equal to the nu- 
merator. The resonance condition as given by Eq. l60l 
implies 



k2d = 2mr; — >■ /c2 = npd 



(61) 
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where n is an integer and the characteristic momen- 
tum scale pd = l-K I d has been introduced. From 
Eq. I55bl and Eq. [HI] we obtain the following condi- 
tion for complete transmission of an incident wave: 



\r? sin^ a + n"' cos^ a( -f 



(V-E)d 



(62) 



or equivalcntly in terms of "renormalizcd" constants 



n[l 



n2(^)2] 

Pa 



PdVa 



(63) 



Eq. 123] gives the resonance condition, either for 
resonant energies En{a,d,V) or for orientations 
a„(d, V — E), for full transmission. 

The limiting cases are a — )- and a — > 7r/2. The 
latter limit corresponds to normal incidence of a par- 
ticle with 'relativistic' Dirac-Weyl dispersion which 
is treated in Ref. [l^ , where it was shown that there 
is complete transmission even if the potential barrier 
is large. The resonance condition for this limiting 
case can be obtained setting a = it/ 2 in Eq. [551 
The a = limit becomes the case of conventional 
massive particle tunneling, whihc must be treated 
separately (see the following subsection). The semi- 
Dirac system provides for, and interpolates between 
smoothly, the two very different limits. FigureUpro- 
vides a schematic illustration where there is a single 
resonant orientation of the barrier. 



C. Limiting case a = 

This case corresponds to the potential being per- 
pendicular to X (the non- relativistic direction), so 
ky~0. The Hamiltonian admits evanescent as well 
as propagating wave solutions only in this case; in 
a sense the relativistic character dominates the be- 
havior except at a^O. It is instructive to follow the 
mixing of the positive and negative energy compo- 
nents. Operating on propagating waves e^*'^"^^ the 
Hamiltonian in Eq. |49] takes the following form in 
the k space: 



2m 



(64) 



same:(£' = 2^). The momenta in regions {/, ///} 
and // are denoted by fc" 



■\/2m|T^^^^ respectively, 
function in the three regions are 



= \/2mE and k!^ = 
The form of the wave 



(65) 



-00 < X < ~d, 



—d < X < d, 



d < X < 00, 

where r",t"',t'(, r'{,t'(',r'{',t 2, r'2 are constants. In 
Eq. for regions / and /// the evanescent waves 
are constructed in such a way that they don't diverge 
when I a; I becomes large. There is no backward trav- 
eling wave in region III. jijp is the transmission 
coefficient. 

Equating the wave function and its derivative at 
the boundaries x = and x = d, for the transmission 
coefficient we obtain 



\t'^\' = \- 



ikn d 



e-^--'^{k'i + ik'lY - e^'2'^{k'^ - ik'{f 



'•(66) 



Eq. [66l is the same as that given by Katsnelson et 
al^ in the context of the tunneling probability for 
the bilaycr graphene dispersion, fcj gets large as 
the potential V gets large. Because of the presence 
of the exponential factor e^^^ <i in the denominator, 
the transmission coefficient given by Eq. IBHl goes to 
zero as the potential goes to infinity. Thus there is 
no perfect transmission when the potential is in the 
non relativistic direction and the particle is incident 
normally, as mentioned above. 



with the conventional massive eigenvalues ±tt^. For 
evanescent waves the eigenvalues are interchanged, 
resulting in a mixing of positive and negative energy 
functions in a way that does not occur with non-zero 
ky. 

The energy of the incident particle for both 
the propagating and the evanescent cases are the 



IX. SUMMARY 

In this paper several low energy properties of the 
semi-Dirac, semi-Weyl degenerate semimetal have 
been studied. Whereas some of the properties are in- 
termediate between the conventional parabolic and 
the linear "Dirac" dispersion, as is the case for the 
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NON-RELATIVISTIC 
DIRECTION 



FIG. 4: Complete transmission for various orientations 
of the potential 



cyclotron frequency, some other properties can be 
distinct and rather unusual. The results for the 
Klein scattering for the semi-Dirac dispersion have 
been obtained for normal incidence on an arbitrarily 
oriented barrier in the 2D plane, revealing that an 
electron can tunnel through the barrier with prob- 
ability one, subject to a resonance condition being 
met, except for the direction where linear dispersion 
does not enter the problem. The extreme anisotropy 
of the plasmon frequency is a distinctive feature of 
a semi-Dirac system. Intriguing behavior for the 
Faraday rotation. Hall coefficient and heat capac- 
ity have been provided. Finally, we note that the 
behavior of the orbital susceptibility is distinct from 
both quadratic and linear systems, being strongly 
dependent on doping level. 
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